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In the framework of zeta-function approach the Casimir energy for three simple model system: 
single delta potential, step function potential and three delta potentials is analyzed. It is shown 
that the energy contains contributions which are peculiar to the potentials. It is suggested to 
renormalize the energy using the condition that the energy of infinitely separated potentials is zero 
which corresponds to subtraction all terms of asymptotic expansion of zeta-function. The energy 
obtained in this way obeys all physically reasonable conditions. It is finite in the Dirichlet limit 
and it may be attractive or repulsive depending on the strength of potential. The effective action 
is calculated and it is shown that the surface contribution appears. The renormalization of the 
1/^ , effective action is discussed. 
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I. INTRODUCTION 
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I Recently there was great interest to Casimir effect for semi-transparent boundaries |3, S ITsl ITol l20l 123^ . It was 
suggested to consider the potential with zero support instead of the rigid Dirichlet condition. This potential models 
the semi-transparent boundary condition. The space is not divided into separate parts by Dirichlet boundaries and 
some modes of field may go through the boundary. Firstly, this kind of calculations were made in Refs. 0, 0, |^ . 
It was claimed that in the limit of infinite strength of potential, which was called Dirichlet limit p^ . the energy is 
, divergent and it does not coincide with that obtained for Dirichlet boundary condition. The authors of noted 
T— I ■ that in framework of QED it is impossible to obtain finite result. This divergence appears in the energy only. The 

' force, which is the derivative of energy with respect the position of plate, is finite '20] in the Dirichlet limit. 
1^ . The same problem was emphasized in Ref. 6]. Singular potential brings additional surface contributions to heat 
kernels coefficients Q and to the effective action and there is no universal method to make unambiguous renormaliza- 
tion procedure and fix all parameters. Necessity to consideration of the surface contributions to energy in framework 
of bag model was noted in Ref. 1 , l§j . From the point of view of field theory the delta- like potential may be considered 
Qh! as brane theory and the parameter of potential strength is regarded as brane's mass In framework of quantum 
^ • field theory this parameter has to be renormalized, too. Renormalization group equation for brane with co-dimension 
greater then unit was obtained in Ref. |l2j | . It was noted that for renormalization we have to introduce into theory not 
^ 1 merely brane mass but also brane tension and others parameters 12] the number of which depends on the bulk 

' action. There is another analogy of system under consideration with field theory in curved space-time with singular 
■ scalar curvature. For example, in space-time of short-throat wormhole the space-time is everywhere flat except 
d ' the throat where the scalar curvature is singular. The strength of the delta-potential corresponds to non-conformal 
coupling constant which has to be renormalized too |22j ] . 

In this paper we reexamine this problem in framework of zeta-regularization approach considering in details some 
simple model system in Sec. ^ namely, single delta potential, step function potential and three delta potentials. 
The main features of these models are summarized in Sec. IIIII It was noted that the singularities appeared in the 
limit of infinite strength are connected with brain (potential) itself. There we suggest a method to extract physically 
reasonable expression for energy which obey all conditions we need. It is shown that the same result may be obtained 
by Lukosz renormalization procedure. The Casimir energy may possesses the maximum and the Casimir force may 
be attractive as well as repulsive. Close analogy with " surface energy" is noted. The effective action in framework of 
zeta-regularization for theory with delta potential is considered in Sec. IIVI Recently similar approach was used 
by Toms in Ref. for free scalar field. The same ambiguous of renormalization procedure in surface term is 
recovered. We observe that this ambiguous is connected with some surface quantity peculiar to the potential (brane) 
itself. By using approach of Sec. IIIII we obtain physically reasonable result which is coincides with that obtained in 
Sec. IIIII We note that the Lukosz renormalization procedure takes off all singularities of effective action. 
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II. THE ZETA-REGULARIZATION APPROACH 

Let us consider the massive scalar field in + 1 dimensional Minkowski space-time with potential V which depends 
on single coordinate x, only. In the framework of zeta-regularization approach 0j ^^e regularized energy is 
defined by following expression 

(1) 



where 



is the zeta-function of operator 



2s 

(«) 

(") 



Dn = -AN + m^ + V 



with eigen- values w^^^ . Parameter with dimension of mass has introduced to keep right dimension of the energy. In 
accordance with Ref. |2| one defines the Casimir energy by relation 

Ei^^^lun(E^^Hs)-EZHs)), (2) 

s^O 

where 



«+i_„r(s + 



mJ 2(47r)^/2 ^ ^ r(s - i) 

is the divergent part of energy, where are the heat kernel coefficients of operator D^. This expression obeys to 
the physical reasonable condition 

lim Ef^ = 0. (3) 

m — ^oo 

The spectrum of this operator is numerated by one discrete number n, and continuous numbers ^n-i G (— oo, +00) 

By integrating the definition of zeta function over continuous variables kjv_i we obtain the relation 

2 (47r)— r(s- 2) 2 

by using which we have to solve now the one dimensional problem only. We note that the using above relation in Eq. 
([T|) gives the energy per unit square of plate, that is surface density of energy. To calculate the zeta-function wc will 
use approach in framework of which the zeta-function is represented in the form below 

C(s,Di) = ^^^^^^ r dk{e ~m')~'^\n^!{ik,R,R') (5) 
TT 7„ dk 

where the spectrum of energy k^ = m? — lo^ is found from relation ^'(fc, R, R') = 0. Let us proceed to consideration 
of different forms of potential. 

A. Single singular potential V = \s5[x) 

We consider a singular potential V — Xs5{x) and two iV — 1 dimensional Dirichlet plates at points x = R and 
X = —R'. Because of Eq. (O we need to consider the imaginary energies only. Let us denote ~ uP' = k"^ and 
consider the one dimensional problem in imaginary axis k — > ik. In this case the eigen-problem equation reads 

+ e + K5{x) </. = o. 
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In two domains x G (0, R) and x £ (— -R', 0) we obtain the general solution of this equation 

(t)+ ^ C+ e^^' + C+ e-^"" , 



Four constants and obey to four homogenous equations 

(/)+(o) = r(o), 

<^+'(0) = r'(O) + A,0+(O), 
= 0, 

(t>^i-R') = 0. 

The solution of this system exists if and only if the determinant of the system equals to zero. This determinant is the 
function 'i>{ik, R, R') which is used in Eq. 

*(zfc, i?, R') = I \e-''^'A' + c'^^'a] , (6) 

k I i 

where 

1 



1 



A = j- [A,e-^-^ -(2fc + A,)e'=«] . 

It is easy to see from Eq. ^ by equating 'i'{ik, R, R') to zero that the boundary states appear if 

, R + R' ,_, 

This boundary state is localized in the potential (brane) and in the case of whole space i? = i?' — > (X) it has the form 
below iia]: 

Because the possible boundary states have already incorporated in Eq. (|3J| we assume arbitrary sign of As . Nevertheless 
we suppose that the inequality 

Xs > -2m (8) 

is valid. In opposite case we can not use the theory in present form (see Ref. 0|). 

To find the heat kernel coefhcients which are defined as asymptotic expansion over mass, 

11°° 1 

^ n— 

we take the asymptotic expansion of zeta function (|42|) over k oo. One has 

1 1„ ij, «')" . + - 1 - i;(5JTX) ^ « + - ^ |(- 5FfT' <") 

By using this expression we obtain heat kernel coefhcients 

Bn = R -\- R' Bi — —\/tt, 

2 

\ \ 2n , — 

As \ a/tt 
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Then, using Eqs. I@J and (0) we arrive at the following expression for Casimir energy 

1 1 

2 47r 2 i ^ + i ) Jm 



where 



r) 1 ^ A' 

ZN[k,Xs,R,R'] = 7^1n«'(fc,i?,i?')-(^ + ^') + T-E(-l)'TOT- 



1=1 



Let us consider some limiting cases. In massless limit m the energy is divergent 



(-1) 



N 



-Bn+1 In 



in agreement with Ref. 4]. The energy in whole space {R ~ R' —f oo) 

1 1 



(11) 



(12) 



(13) 



S(^'[A^,i?-> oo,i?' ^ oo] 



2(47r)f r(f + 1) 



dk{e - m^f 



A. 



N 



(-l)^+im^ / A 



4(47r)^/2 V2to/ lr[^] 



N+l 



iFi 



I N + i ( \s 
2m 



^' 2' 2 



k{2k + \s 

Xs 1 



EM) 



(=1 



2'fc'+i 



(14) 



2m r[^] 



1,1: 



iV + 4 / A 



2 ' V 2 



is finite but it is ill defined in the Dirichlet limit, A^ — s- oo. Here the 2F1 is the hypergeometric function. The leading 
divergent term in this limit coincides with that in Eq. (|13|l . The energy is well defined in the limit A^ — s- 0. In this 
case 



ZN[k,\s^O,R,R'] = 



2{R + R') 

,2k(R+R') _ I 



which corresponds to the Casimir energy for two Dirichlet plates at points —R' and R. 

Let us consider the energy for particular dimensions A'^ = 1 , 3 in manifest form. Then for A'^ = 1 we obtain 



dx(x2-l)i/2Zi[mx,A„i?, i?'], 



where 



Zi[k,Xs,R, R'] 



= ^^8k^{R + R') + 2kXs 
X {2fc2 [2fc(e2'=(^+«') - 
In the limit R = R' ^ 00 



-2 + e^'^^ + 6^'=^' + 2fci?'(e2'=« _ 1) + 2kR{e'^''' - 1)] + (e^''^^ - l){e'^^' - 1)A2} 



l) + (e^^^^'-l)(e 



2kR 



L)A.]}" 



£;(i'[A,,i? ^ cx3,i?' -> 00] 



47r 



As /, A2 A, 

TT \j A arccos - — 

2m 



This expression coincides exactly with that obtained in Ref. [13 ■ It is ill defined in the Dirichlet limit A^ — > 00: 

i?W[A, ^ 00, i? ^ 00, i?' ^ 00] « ^(1 - In^) - ^. 

47r 2m 4 

The force acting on the plate coincides with that obtained in Rcf. 20] in massless case. We have to take limit 
R' ^ oo and then take the derivative with respect R with sign minus. After integrating by parts we obtain the force 
in the form below 



(15) 
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To compare with Ref. [23 we have to take Umit A' ^ 00 in Eq. (2.13) of this paper and change the notations of 
parameters a ^ R and A Agi?. As was noted in Ref. this expression is well defined in the Dirichlet limit 
Xs 00. 

In three dimensional case, iV = 3, one has 



£;f)[A3,i?,i?'] = --^ / da;(x2-l)3/2z3[mx,A„i?,i?'], 



(16) 



where 



Z3[k,Xs,R, R'] 
^^32k^{R + R') + 8k^X, 

2kXl 



-2 + e^'^^ + e"^^' + 2A:i?'(e2'=« - 1) + 2A:i?(e2'=«' - 1 



-2 + e^''^^' + e^^""' 



+ (e^'^^ - l)(e2^-^' - 1)A^} {8fc4 [2fc( 



Ak^Xl 

2k{R+R') 1^ I /„2fc7?, i^/„2fe_R 



-2 + e^''^ + e^'^^' 



l) + (e^««_l)(e^««_l) 



In the limit R = R' ^ 00 we obtain the following expression: 



E^^^[Xs,R-^ oo,R' 00] 



A. 



Ac. „ A „ 



5767r2 

which is divergent in the Dirichlet limit As — > 00 



247r - 24— - 97r^ + 8^ - 6(4 arccos 



A, 

TO 



) [As ^ 00, i? ^ 00, i?' ^ 00] 



A„ TO A: TO^A, 



1 



727r2 647r 327r2 247r 167r2 



4 + "^'0ln-. 
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B. A single step function potential 

Let us regularize delta function by step function by relation 

J(.) = lim(0' N>^' 

and consider the energy for finite value of the regularization parameter e. In this case we get the following expression 
for function for R, R' > e: 



^°"'{ik,R,R') = - 



where 



A[ = e'=«sinh(2e/fc,):^^^-e-'=^+2^'^' 



sinh(2efce) 



kR-2kt 



p + fc2 

sinh(2eA:<:)-^— h 2cosh(2efce 



fc2 + fc2 

-,-kR 



e-'"''smh{2ek,) 



2 cosh(2efce 



and k1 = k"^ + Xs/2e. In the limit e this function tends to that obtained in last section and given by Eq. ©: 

lim ^°''\ik, R, R') = ^{ik, R, R'). (17) 

C-+0 

The first three heat kernel coefficients are the same as for delta potential case, the rest coefficients are divergent in 
the limit e ^ 0: 



K+i = (-1) 



2r(n+ i) 



(18) 



7,1 



^ I e-'\ n=l,2,3. 
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To find the function ^ inside the potential we have to assume R,R'<e and consider additional Dirichlet boundaries 
far from the potential at points x = ±H. In the end we tend these boundaries to infinity: H ^ R, R, e. Because we 
have three domains confined by Dirichlet boundaries the function 'i' is the product of three functions: 



'i>T{tk,R,R') = ■fl^^^'if^ 



(19) 



where 



sinh(fce(i? + -R')) 



*3 = ^ R)_ 

The heat kernel coefficients have the form below 



Bl 



B 



n+i 



2(-l)" (\ 



(-1)" 



-n+l 



(20) 



2r(n+i) 



y ) e-", n= 1,2,3,. 



The difference with above calculated heat kernel coefficients appears at terms with half-integer indices: 

- 5°:*. - (-l)"+^2V^ {^^y e-\ n = 0, 1, 2, . . . . 

In the framework of zeta-regularization approach the energy has the following form: 

1 1 Z""" 

i?W[A.,i?,i?',e]=--— - / dk(k^-m'fl^Z^{k,\,,R,R!,^\, 

2 47r 2 1 + 1) Jm 



where 



Z°j^\k,K,R,R! 



dk 



\n^°"\k,R, R') - {R + R') 



- 



1=1 



Z';^[k,K,R,R',e] = 



(-i)'r(; + i) Aj 

3 



1 

' k 
[f] 

E 

1=2 



(21) 



(22) 



(-l)'A 



,£(i+i^ 



- 1 



(23) 



(-i)'r(z + i) Ai 



i\ 



E 



(~l)'Ai 



, £([+!) 



In the last expression the limit _ff — s- oo is assumed. Because of relation (|17|l this expression and the energy are 
divergent in the limit e ^ starting with dimension = 3. In the limit of whole space R = R' = L/2 oo, the 
energy is finite. It is ill defined in the Dirichlet limit, A^ oo, or/and m — > 0: 



N 



(-1) o , 

——w+tBn+i In — 
(47r)^~ ^ As 



(24) 



where the heat kernel coefficient is given by Eqs. (|18|l and H20I) . 

The energy calculated for this model potential has additional peculiarity for R, R' w e. The energy is divergent in 
the limit i? — > e or i?' ^ e, but it is finite for R ^ e or R' ~ e. Indeed, let us consider the asymptotic expansion of 
function 
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The half-integer coefficients obtained by this expression are different from (|18() starting from | coefficient: 

Bo = 2e, Bi = 



Br, 



2(-l)" (\ 



(25) 



B„ 



(-1) 



ri+l 



71 = 1,2,3, . 



It is easy understand the origin of this problem. To obtain the heat kernel coefficients H18|l we assumed that R,R'>e 
and for this reason we turn down all terms as e~'^^^~'^^ and e"^^^ and after this we set R = R' = e. We may turn 
down these kind of terms for arbitrary small, but non-zero, value of i? — e and R' — e. 
expression for asymptotic expansion 



We obtained the following 



dAn^':"'izk,R,R%:B'=, 



2 
"k 



2 



k 



2ke 

"fcT' 



If we put R — R' — e at the beginning we obtain another expression for asymptotic expansion 

^.An^ri^k,e,er = -^+^-^. 



The difference is not zero 



d,\n^°"'{tk,e,er = -- 
k 



2 

kl 



2e/c3 



To reveal this divergence in manifest form let us use another form of asymptotic expansion. We keep all terms as 
g-fc(ij-e) g^jj^ g-fc(-R -e) g^j^jj consider asymptotic expansion with these terms as constants: 



^« r -2fe(_R.-c) 



{R + R') 

2k{R'-e)' 



A. 



4efc2 



2e- 



-2k{R- 



"^R - e) - e-^^^^' -^\r' - €) 



32e2A;4 



From this expression we observe that additional terms appears. The contributions of them to even degree of k (heat 
kernel coefficients with integer indices) are insufficient. For R,R'^e they are exponentially small, but for R,R' — e 
they are zero. The contributions to odd degrees of k (heat kernel coefficients with half-integer indices) are important. 
They are exponentially small for R, R' ^ e, but they are constant for R, R' = e. The first non-trivial contribution 
starts from l/fc^. Therefore this kind of expansion reproduces right the expansion for R, R' e as well as for R, R' = e. 

Let us use this expansion to calculate the zeta-function and the energy. It allows us to find in manifest form the 
divergence as the boundary position R close to boundary of potential e. In the 3-dimensional case the divergent 
contribution is due to the following term 



A. 



127r2e 



(fc2_rn2)3/2 



dk. 



For i? — e <C 1 we obtain the following behavior of energy close to boundary of potential 

A.. 



127r2e|i?-e| 



(26) 



The application of this approach to internal function l|19|) gives the same result but with opposite sign. 

It is obvious that this divergence connects with model under consideration. If we adopt another kind of regularization 
for delta-function as a consequence of analytic functions we obtain no divergence except for infinitely close position of 
plates each to other as should be the case for Casimir force. For position of Dirichlet plates exactly on the boundary 
we obtain the following finite expression for energy 



1 r°° 
127r^ 



-2e 



k 



2fce w , ^ As 

— COth(26fc,) + — 



A. 



3A2 



2efc3 16efc4 



(27) 



C. Three singular potentials V = Xs25{x + R') + \sS{x) + Xsi5{x — R) 

Let us consider three delta-potentials 

V{x) = Xs26ix + R') + \s5{x) + \si5{x - R) 

and the field with Dirichlet boundary condition at surfaces x — —I, L, where I, L > i?', R. 
The ^ function has the form below 

k(l+L) , 

i?, R!, L, Z, A,, A,i, A,2) - — ^ {- (2fc - A,) (2fc - A,i) (2fc - A,2) e-2'=('+^) 

- (2fc - A,) A,i (2fc - A,2)e-2fe('+^) + A,A,i (2fc - A,2)e~'"('+^-^^ - A, (2fc + A,i) (2fc - A,2)e-2'=' 

- (2fc - A,) (2fc - A.i) A.2e-^K^+^') + A. (2fc - A,i) A.2e-^K'+^-«') + A.A.i A,2e-2H^+«'-«) 
+ A.A,iA,2e"''('"'''+'') - (2fc - A,) A,iA,2e"'''('''+'') + (2fc + A,) A,iA,2e"'"('+''"'''"'') - A, (2fc + A,i) A. 
^ (2fc + A,) (2fc + A,i) A,2e-''-('-^') - A, (2fc - A,i) (2fc + A,2)e-2'-^ - A,A,i (2fc + A,2)e^2fcfl 

^ (2fc + A,) A,i (2fc + A,2)e-2Mi-i?) + (2fc + A,) (2fc + A,i) (2fc + A,2)} • 

For — = this function coincides to that considered in the first section. The following relation is valid: 

*(fc, i?, i?', L, 1, 2X„ A,i = 0, As2 = 0) = ^'(fc, i? = 0, i?' = 0, L, A^ = 0, A^, A^) 

as should be the case. But the energy does not obey this kind of relation. Indeed, let us consider the case of two 
potentials and put A^ = 0. One has 



26 



-2kR' 



*(fc, R, R', L, I, Xs = 0, A,i, A,2) = {-{2k - A,i) {2k - A,2)e-'"'^ - A,i {2k - A,2)e-2'^(^+^) 

- {2k Xs,) A,2e-^'=(^+^') - A,iA,2e-^'=(^'+^) + A,iA.2e-^'=(^^-^'-^) (28) 

- {2k + Xsi) A,2e""'(^"'^') - A,i {2k + A,2)e-'^'(^-^^ + {2k + A,i) {2k + A,2)} • 

In this case 



2kL 



l^sl + As2 



ZN[k,0,Xsi,Xs2,R,R'] = ^ln^{k,R,R',L,l,0,Xsi,Xs2)-{L + l) + ^-J2{-l) ^z^^+i 

n 1 ^ A' 

ZN[k,Xs,0,0,R,R'] = —ln^{k,R,R',L,l,Xs,0,0)-{L + l) + --Y,{-^y^r^- 

1=1 



The difference has the following form: 

N 



AZ^-J2{-1) 



2X1 - (2As 



1=2 



and it is non-zero starting from dimension iV = 3. This fact originates from renormalization procedure. 
In the limit of whole space, L — I oo, we get from (|28ll 



2k L 



^'(/c, R, R', L, I, Xs = 0, A,i, A,2) = {-A«iA,2e-''=(^'+^) + {2k + A,i) {2k + A,2)} ■ (29) 

If we then put potentials to infinity R ~ R' ^ oo we obtain non-zero result 



N 



Xsi Xs2 sr^, iwAgi+A' 



Z.ik,0,X.,,X.2,R,R'] = -^(^ri^^ - ^(^ri^^ 2',..: 

\i Ail + Ai 



s2 



- - E (-1)' 



^s2 



l=N+l 
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This is in contradiction with expected result that the energy has to be zero. We observe that in this Hmit the energy 
is the sum of energy of two single plates. Therefore we conclude that this energy is connected with potential itself and 
it has the additivity property provided the infinite distances between potentials (without interaction). For J plates 
in the limit of infinite distances we get 



2- = -E' 



A. 



J N 



A 



1=1 ^ 1=1 1=1 

Therefore, with each j-th singular potential \sj5{x — Rj) we may connect the "surface" energy by relation 



2(4^)f r(f + 1) A 

This energy is ill defined in the limits m ^ or/and Xsj 



\N/2 



A. 



N 



mx{2mx + Xsj) 



1=1 



(30) 



III. RENORMALIZATION PROCEDURE 



Let us summarize useful information from above section. First of all, we should like to note that the divergence 
of above energies is connected with renormalization procedure only. The regularized expression for energy Q , Q is 
well defined in the limit A^ ^ oo or e ^ 0. The renormalization procedure brings a problem. The point is that for 
renormalization we subtract finite number of terms (see (|12|l or 1)22(1 ') of series expansion over small value As/fc ^ 1 
or As/efc^ ^ 1. After renormalization we use expressions obtained for As — > cxd or e ^ 0. Obviously that the main 
divergence originate from the last term of truncated series. Indeed, let us consider the asymptotic expansion over k 
the integrand 

i in K. «r- ^ + «■ 4 - + 0(e-, (31) 

The last term is finite in the limit As ^ oo and the difference 

— ln-^ik,R,R')-—ln-^ik,R,RT (32) 



is finite in this limit, too. But if we use truncated series expansion 

1 1 ^ 

-ln*(fc,i?,i?')- = i? + i^'-^+EM) 2^,H 

1=1 



the difference H32f) is infinite in the Dirichlet or massless limit. The main divergence comes from the last term which 
gives the contribution to heat kernel coefficient B n+i which in turn gives the conformal anomaly. 

Second, with each potential connects the inner quantity given by Eq. l|5n|) which we called "surface" energy. It is 
defined as energy of single delta-potential provided the boundaries are in infinity. This quantity possesses an additive 
property: the " surface" energy of set delta-potentials is the sum of " surface" energies of each delta-potential provided 
the infinite distance between potentials. We note also that this quantity is connected with potential itself. In self- 
consistent theory which takes into account gravitational field of brane |25( , this energy will give some contribution to 
gravitational field of brane. 

For this reason let us consider another renormalization procedure. Instead of extracting terms which will survive in 
the limit m — + cx3 we subtract all terms in asymptotic expansion of zeta function. It is easy to see that this procedure 
corresponds to subtracting "surface" energy of delta potentials, which is defined for R, R' oo. It seems reasonable 
that the Casimir energy is zero if all potentials are in infinity. Indeed, all terms of asymptotic expansion in H31|l may 
be obtained by taking the limits R, R' ~* oo: 

I m *(^, R, iru^.^^ = ^ + ^' 4 - ki2txr) 

It is easy to see that the condition Q is still valid for this kind of renormalization. In this case the energy for single 
delta-potential considered in Sec. Ill Al has the form Hll() 

i?W[A.,i?,i?']=-— ^— - / dk{k'-mX/^Zr,[k,Xs,R,R'], (33) 
2(47r) 2 i (^ + 1) J„ 
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where function 

2 f 4fc2(i? + R') + X, \-2 + e2'=-R(l + 2kR') + e^^^' (1 + 2kR) \ + \{e^^R' -\)R+ (e^^^ - l)R' 
ZN[k,Xs,R,R']^ ^ ^ ^ ^ 



{2k + Xs) [2k{e'^HR+R') _ l) + (e^kR _ i)(e2fci? _ i)x^] 

(34) 

does not depend on the dimension of space. It is easy to see that the difference of these energies is that divergent 
expression (|14|l which we called "surface" energy H30() : 

[Xs,R, R'] - Ei""^ [A„ i?, R'] = Ef^ [A„ i? -> oo, i?' -> oo] = Ef\ 

The renormalization procedure suggested by Lukosz [l^ (see also 8]) gives the same result. Indeed, according with 
this approach we have to include our system to large box, for example with x — ±_ff, where H > R, R' . The space is 
divided into three domains: Mi : {a; G [—H,—R']}, M.2 ■ {x G [—R',R]} and M3 : {x e [R,H]}. The renormalized 
energy is given by equation 

= E(f\Mi) + Ei^HM2) + E^fHMs) - E^f^M), 

where the last term is calculated for whole space but without internal boundaries. In the end we have to tend H —^ 00. 
It is very easy to apply this formula in our case. Because in domains Aii and AI2 there is no potential we may use 
our formulas with As — 0. Therefore we obtain 

Ai?(^) = Ei^^ [A, = 0, -R', H] + Ei^^ [Xs,R, R'] + sf^ [A, = 0, H, -R] - E^^^ [X,,H, H] 

_ 1 1 r, .,2_ 2.3/2 fc^vl'(»fc, ~R', H)^^=o^{tk, R, R')^{ik, H, -R)x^=o 

" 2(4^)f r(f + 1)X„ '^^'^ dk'"" ^{^k,H,H) 

In the limit — > 00 we obtain the expression for the energy which coincides exactly with that obtained above and 
given by Eq. ((23 • 

Let us consider some special cases of energy obtained. In massless limit the energy is finite for any dimensions. In 
the Dirichlet limit As — > 00, 

'2R 2R' 

Zpflk, Xg —* 00, R, R'] = 



a2kR _ I ' g2fei?/ _ I ' 

the energy, as expected, is finite and it is the sum of two Casimir energies for domains (— i?', 0) and (0, R): 

^(^)[As ^ 00, R,R']= )[As, 0, R'] + Ef^ [As, i?, 0], 
as should be the case. For As = 

2{R + R') 



Zjv[fc,As =0,i?,i?'l 



a2k(R+R') _ I 



which gives the Casimir energy for two plates in points —R' and R. Therefore, this energy has satisfactory behavior 
for all reasonable cases. In one-dimensional case N — 1 the force acting for one plate (i?' 00) due to singular 
potential coincides with that obtained by Milton in Ref. in massless case and it is given by Eq. (|15|) . This fact 
is evident because the difference of two energies is constant which does not depend on the distance to boundary. 

The dependence of the Casimir energy H33() for three dimensional case is plotted in Fig. ^for i? = i?' = L/2 as a 
function of Lm for different values of As and the fixed mass m. We observe that the energy is negative for As > 0. It 
is easy to see this statement from Eq. (|34f) : in this case all terms of Z are positive. The slope of curves is negative 
for all position of plates. It means that two plates are attracted as in usual Casimir case. For negative As a maximum 
appears at position Lm and the energy in the maximum growths with increasing |As|. For distances smaller then L^ri 
the plates are attracted but for L > Lm the force is repulsive. The repulsive character of Casimir force due to delta 
potential was observed in Ref. [l^. Because of the present theory is valid for Xg/m > — 2 we obtain from Eq. O that 
for Lm > 2 the boundary state appears. The energy of this state E = \Jm^ — fc^ is found from following equation 

kL 2k 
tanh — = — — . 
2 Xs 

It is easy to show from Eq. (|34|) that the energy is positive for large distance between plates vnL 00 and for 
> Xs/m > —2. Therefore for large distances between plates the Casimir force is repulsive and boundary state 
appears localized at the position of potential. There is a close analogy of this point with widely discussed now 
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FIG. 1: The Casimir energy for two plates in positions x — ±L with singular potential V = \s5{x) between them for As = —1.9m 
(thick curve), As = (middle curve), and As — 4m (thin curve). 



"surface states" [Tj|. It was shown the repulsive character of "surface states" which may be regarded as boundary 
states Q. 

In the case of step function potential the Lukosz regularization or renormalization procedure suggested above give 
the following expression for energy {R, R' > e) 



£;(3)(A„i?,i?',e) 



127r2 



dk(k^ - m2)3/2__in 
ok 



d k'^^l^'\ik,~R',H)x,=o-^T\ik,R,R')^f'\ik,H,~R)x,^Q 



(35) 

For the case of position of plates inside the potential, R,R'<e we put the system in large box with size 2H. We 
have three domains: x e {—H,—R'), x G {—R',R), x G {R,H). For renormalization we have to subtract the energy 
without internal boundaries at R and —R'. This energy is defined by function vpo"* because the boundaries at ±H 
are outside the potential. Therefore one has: 



E^^\Xs,R, R',e) 



dk{k^ 



,2^3/2 



d ^^ k'^-^f{ik,R,R') 
dk ^°"*{ik,H,H) 



(36) 



Due to relation H17|) it is obvious that in the limit e 
calculated above for S potential 



In the limit A, 



— !■ oo one has 



the energy ()35|l is finite and it coincides with energy 
^(3) (A,, i?, i?', e ^ 0) = Ei^\Xs,R, R'). 

2{R-e) 2(i?'-e) 



127r2 



dk{k^ ~ m^) 



2\3/2 



g2fc(i?,-£) _ I g2fe(J?/-e) „ I 



(37) 



which is exact the sum of two Casimir energies for domains x e [—R', — e] and x € [e, R] as should be the case. For 
zero potential As = the expression 1)35(1 reproduce the Casimir energy for two boundaries at points x — ~R' , R. 

The numerical calculations of Casimir energy is reproduced in Fig. [21for R = R' = L/2, As = 1, e = 0.1. We observe 
the divergence of energy at the boundary of potential L = 2e = 0.2 which was noted in Sec. IIIBI For all position of 
Dirichlet plates the Casimir force is attractive. In the right figure we show details of the energy close to this boundary. 
The smaller e the closer the position of this divergence to origin and the closer external part to the case calculated 
for delta- function case. For small e and negative As the Fig. ^is reproduced. 

In the case of three delta functions in framework of this procedure we obtain the following expression for the function 

Z: 



^{k,R, R',L,l,Xs,Xsi,Xs2) 



(2A: + As)(2fc + Asi)(2/c + As2) 



The same expression may be obtained by Lukosz renormalization procedure. We have to use this approach for each 
solitary potential provided that others are in the infinity. This expression is well defined in all limits we need. For 
example, in Dirichlet limit As,Asi,As2 oo the energy is the sum of fore Casimir energies for plates in points 
X = -l,-R',0,R,L. 
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FIG. 2: The Casimir energy for 3-dimensional case and for R = R' — L/2,Xs — l,e = 0.1. The boundary of potential is at 
point Lm = 0.2. There is a divergence in this point which is plotted in detail in the right figure. 



IV. EFFECTIVE ACTION POINT OF VIEW 



Let's consider the same problem from the point of view of effective action. One has the scalar field ip living in 
(3 + 1) dimensional space-time with singular potential concentrated at the surface E. To make clear regularization 
procedure we consider a nonlinear theory with the following action 

^^Um (^^^^' ^ ^ + ^^^'^ ^^^^ 

Here the A is cosmological constant and A^ is surface tension. As was noted in Refs. j^, we have to include the 
surface terms at the beginning for correct renormalization procedure. In fact, the action (|38|l describes a brane with 
co-dimension one in four dimensional space-time with tension As and mass of brane . The renormalization group 
of brane with co-dimension more then unit was considered in Ref. 12]. Similar problem for free scalar field without 
self-action was considered in Ref. |24| . 

To calculate the energy of quantum fluctuation we consider the field at finite temperature T = 1/(3 and in the end 
of calculation we take the limit of zero temperature /3 oo. For this reason we proceed to the Euclidean regime 
t —f —IT with T e [0, /?]. Let us consider the Casimir problem - two Dirichlet plates at points x = R and x = —R'. 
The singular potential is concentrated at point a:; = 0. 

We divide the field on classical part and quantum fluctuations: ip — ipd + In one loop approximation and 
framework of zeta regularization approach we have the following expression for action 

where 



D4 = -A4 + + Xipli + XsS{x). 



(39) 
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In the limit s — s- 



where /i^ = /i^e ^. 
Due to the relation 



C(s,i?4) = 



1 C{s-1,D2) 



An 



s- 1 



(40) 



we consider two dimensional eigenvalue problem 

dx dr 



(j) = UJ (j) 



with Af^ = + '^s^ci- Taking into account the periodicity over r we obtain 



d f2im 



0, 



where n = 0, ±1,±2, . . .. Let us denote {j-^^ + — cj^ = and consider the solution of above equation in 
imaginary axis p — ik. In this case the equation reads 



-Q^+k^+^sSix) 



= 0. 



(41) 



This equation has already considered in last section. Using the results we arrive at the following expression for zeta 
function 



sin(7rs) 



M y r dk{e-M, 



n— — oo 



''-'^\n^{k,R,R'), 



(42) 



where = + {21171/ (if. 

Taking into account this expression we obtain heat kernel coefficients 



2n-l 

/ A ^ \ 

B 



n= 1,2,3,...., 



where J7 = {R + R')(3 is the volume and S = /3 is the square of single boundary surface. These heat kernel coefficients 
may be obtained by multiplication to S the coefficients obtained for delta potentials (|10|) . We observe that the zeroth 
coefficient is volume of the system and all others coefficients are proportional to square of boundary. To note this 
moment it is better to define the density 6„ of heat kernel coefficients by relations Bq ^ boQ and i?„ = 6„E, n > 0. 
In calculations of heat coefficients the asymptotic relation ^ 



+00 



1 r(^ + |-^) „,i- 

T{s + l) 



was used. 

To obtain the effective action we have to calculate Cl^l? D2) and C'(^li ^2)- For this aim we will use the formula 
(|42|l . We subtract from and add to integrand the first five terms of asymptotic expansion of it and represent the zeta 
function in the form 



sin(7rs) 



(43) 



^ V / dk{k'-M^)-^Z,[k,Xs,R,R'] 



n— — cxD 
4 



47r r(s) ^2 ^ ^2' 

^ ' ri=0 



14 



Using this equation we calculate the zeta function at point —1: 

Ca(-1,^2) = 0, 

and its derivative in this point 

poo 



a(-i,^2) 



n— — oo 

In r 1 



4tt 



1 ( 1 



An 4 



To calculate CLi^^^ ^2) we use the Abel-Plana formula 

^ 1 



f[iy 



3 



2\/nmb3 — 62 



-dy. 



n=l 

In the limit /3 ^ 00 we get 

C(-l,i?2) = -^^ / dx(a;2-l)3/2z3[mx,A„i?,i?']+0(e-"^). 

Then, using Eq. H40|l we express four dimensional zeta function in therms of two dimensional: 

1 



C(0,Z?4) 



C{0,D,) = -j-C{~l,D,)-j-C{-l,D,). 

Taking into account these equations and putting all terms together we obtain the following expression for effective 
action 



s = ns, 



eff 



(44) 



where 



eff 



bo 



ceff 



GAtt^s 1287r2 6Att^ 
M%i - 62 M^bi/2 M%i 



3M4 , M2 

+ „ . o In — 5- 
M 

Mb 



-^2 2 '^4 
+ ^ci + -^'Pcl 



2 ' 



'3/2 



327r2s 247r3/2 327r2 167r3/2 
,2 



■In— + -A, + -A,(p2,, 



327r2 



and i5c is given by Eq. (|16|l . and A-P — ni^ + A(p2j. The last term in Eq. (|44|l can not be regarded as volume or 
surface contribution. 

From above expression for effective action we observe that one has to renormalize not merely parameters of volume 
part of action but also parameters of surface contribution, too. To renormalize effective action we use minimal 
subtraction procedure and make the following shift of constants 



m + 



167r2s 



, A-^ A + 



3A2 
167r2s ' 



A-> A- 



327r2s' 



A., 



A. 



A61 
167r2s 



A. 



A, 



ni^bi — 62 
167r2s 



(45) 



The renormalization of the surface parameters depends on the details of potential (brane) which are encoded in heat 
kernel coefficients. The same renormalization condition was used in Ref. 0]. We may realize our model as a model 
of scalar field in curved space-time but with singular scalar curvature as, for example, in Ref. for wormhole 
space-time. In this case the parameter A^ is proportional to radius of wormhole's throat a and As plays the role of 
non-conformal coupling A^ ^ and it has to be renormalized, too, as was noted in Ref. Therefore the effective 
action takes the following form 



ceff 



3M^ hd"^ 1 2 2 A 4 1 , 

H ^ m — TT H — rn <^^; H '^^-i H — A, 

1287r2 647r2 2 12 ^'^^ 2 ' 

T2" 



,^ff _ APbi/2 A'Pbi M63/2 Ad%i~b2 



247r3/2 327r2 167r3/2 



32^2 
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To fix the values of parameters of model both volume and surface parts we have the only free parameter /i. To remove 
the renormalization unambiguous in effective potential we may use the standard relation for effective potential 

S^Jfm'^m\ (46) 

by using which we obtain that h\m? / ji^ = 1. This value of ji we exploit in surface part of effective action. Therefore 
there is no way to fix surface parameters of the theory and ambiguous in renormalization procedure in surface term 
appears which was noted in Ref. Q . The similar prolalem appears in Ref . ^3 where effective action was calculated 
in space-time of cosmic string. We arrive at the following expression for action for zero value of classical field 



1 


A, 




1287r2 ^ 2 






m63/2 


247r3/2 ' 


167r3/2 + 




A? 


247r ^ 647r 


1927r2 



1 , 

+ 2 A. 

where the last term is defined ambiguously up to finite terms. We may make consequent arbitrary finite renormaliza- 
tion. To avoid this ambiguous let us, in spirit of last section, subtract from effective action all terms which will survive 
in the limit of infinite distance between plates, that is the "surface" energy. It corresponds to reasonable condition 
that the energy in empty space is zero. In this case the energy of the system is 



F- ^ 



5(0)- hm 5(0) 



= sf)[A„i?,i?']. (47) 



We may use another approach to calculate zeta function at point s = — 1 and subtract from integrand all terms of 
asymptotic expansion 



C(s,Z?2) = Us,D2) + Us,D2) (48) 
sin(7rs) 



V / dk{k^ -Ml)-'Z^[k,\s.R.R'] 



n— — oo 
oo 



2 ^ O / 



47rr(s)^ = ' 2 

^ ' n=0 

In this case additional contribution appears as surface term in action: 

S = nS^" + ^S^Jf -f l3Ef^ -f 0(e~"''), (49) 

where 

S'J^iQ) = S^Jf (0) + Ef^ [A,, i? -> C30, i?' ^ oo]. 

The last term is exactly "surface" contribution i^^^'' given by Eq. 1)30(1 . In the limit As ^ oo it is divergent 

?:eff,r.^ ™^ m?\s 5A^ 1, 1 / A^ 9 ^ \ , A^ 

S^'^iQ) « ^ -I ^— + -As H ^ +m^Xs In—. 

^ ^' Utt 327r2 5767r2 2 167r2 I 6 m 



The Lukosz renormalization procedure considered in above sections gives the same finite result without a renor- 
malization. We will apply this procedure to action H44|l before renormalization of constants (|45() . In accordance with 
this approach we have to include our system to large box, for example with x = ±H, where H > R, R'. The space 
is divided into three parts: Mi : {x G {—H,—R')}, M2 ■ {x G {—R',R)} and M3 : {x € {R,H)}. The contribution 
to energy from boundaries is given by equation 

AS = S{Mi) + S{M2) + SiM-i) - S{M), 

where the last term is calculated for whole space but without internal boundaries. In the end we have to tend H —^ 00. 
It is very easy to apply this formula in our case. One has 

AS = 2psr{X. -0)--l^r dxie - ^^r^ljn ' i^)...oy^i^^y ^(fc, H, -R),^-_, 
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where 



tf M^bT /■) 1 1 



In the hmit ^ cx3 we obtain the foUowing expression 

AS = 2(3S^^{\s = 0) + pE'i'^ [Xs,R, R'] 

and Ei^\Xs,R.R'] is given by Eq. (|33|l . Therefore we observe that this procedure takes off a vohime part, aU 
pole terms and the renormalization parameter fi. Therefore, in framework of this renormahzation the energy (with 
dimension surface energy density) has the foUowing form 



E : 



DAS 



5/3 



/2 



127r3/2 



Ks + Ef\K,R, R']. 



First two terms gives tension of brane with quantum corrections. The last term completely coincides with that 
obtained in Sec. Illll bv subtracting all terms of asymptotic expansion and given by Eq. (|34|l . 



V. CONCLUSION 



Let us summarize our results and observations. In framework of zeta-function approach we considered three model 
potentials namely, single singular potential, step function potential and three singular potentials. We note problems 
which were appeared in connection with these potentials: (i) the Dirichlet limit, ^ cxd, is ill defined and the energy 
is divergent in this limit, (ii) the energy calculated for step regularized delta- function is divergent in the sharp limit, 
e ^ 0, when the step function tends to the delta function, (iii) in massless limit to — > the energy is divergent, 
too. All of these divergencies has the same structure - the energy is logarithmical divergent with second heat kernel 
coefficient as a factor (for 3 + 1 dimensional case). 

In the model of three delta potentials we observe that the energy of this system for infinite separation of potentials 
transforms to a sum of energies each of which is peculiar to single potential (brane) itself which we called " surface" 
energy. To obtain physically reasonable result we suggest to define the Casimir energy as energy without this " surface" 
contribution. In this case the Casimir energy defined in this way for infinitely separated potentials (empty space) is 
zero. The same result may be obtained by using Lukosz renormalization procedure. We showed also, that this kind 
of renormalization correspond to subtracting from zeta-function all terms of its asymptotic expansion. The problem 
appears if we truncate this series and subtract the finite terms of series. The Casimir energy calculated in this way is 
well defined for all physical situations and in the Dirichlet limit the delta-function transforms to Dirichlet boundary 
as should be the case. 

Next, we considered in details the effective action for single delta potential in (3 -|- 1) dimensional case in frame- 
work of zeta-regularization approach. To consider renormalization we use the non-linear 0^ model with brane part 
containing the brane tension and brane's mass. The effective action has surface contributions except the volume part 
(effective potential). The renormalization of the surface part is ambiguous. It is possible to subtract pole terms by 
renormalization of the brane parameters. But there is no universal way to fix these parameters. The application 
the Lukosz renormalization procedure takes off all singularities of the model and gives the Casimir energy which is 
coincide with that obtained by suggested procedure. 
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